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Abstract 

Axions are pseudo-scalar particles, those arise because of breaking of Peccei Queen (PQ) symmetry. Axions 
have a tree level coupling to two photons. As a consequence there exists a tree level coupling of axion to 
photon in a magnetic field. However, in an external magnetic field, there exists a new loop induced, axion 
photon vertex, that gives rise to axion photon coupling. The strength of the tree level axion photon coupling 
in magnetic field is known to be model dependent. However in a magnetic field, the new loop induced coupling 
has some interesting features. This note discusses the new axion photon vertex in a magnetized medium and 
the corrections arising from there. The magnitude of the correction to axion photon coupling, because of 
magnetized vacuum and matter is estimated in this note. While making this estimate we note that the form of 
the axion photon vertex is related to the axial polarization tensor. This vertex is shown to satisfy the Ward 
identity. The coupling is shown to have a momentum dependent piece in it. Astrophysical importance of this 

extra modification is also pointed out. 


I. INTRODUCTION 

Axions play an important role in the conceptual aspects of particle physics today. They are believed to 
be associated with spontaneous breaking of global Chiral symmetry U{1)pq (Peccei Queen symmetry), 
postulated to provide an elegant solution to strong CP problem [1,2]. They are ultralight pseudo-scalar 
field [3]. In the Weinberg-Wilczek-Peccei-Queen model ( original ), the symmetry breaking scale was 
assumed to be around weak scale, fw Although the original model, associated with the spontaneous 
breakdown of the global PQ symmetry at the Electro Weak scale (EW) fw, is excluded experimentally, 
modified versions of the same with their associated axions are still of interest; where the symmetry 
breaking scale is assumed to lie between the EW scale and 10^^ GeV. Since the breaking scale of the 
PQ symmetry, fa, is much larger than the electroweak scale fal^fw, the resulting axion turns out to 
be very weakly interacting (coupling constant ~ ff^), very light {rria ^ and is often called “the 
invisible axion model”. 

Till date, it remains elusive to experimental confirmation, however there have been some efforts to 
constrain its parameters through various cosmological or astrophysical considerations. For instance, 
cosmological observational constraints put bounds on its mass (such that the universe is not over 
closed). Through such arguments, the allowed range for the axion mass rria has turned out to be [4-7], 

10"® eV < ma < 10-2eV. (I.l) 

Apart from the one mentioned above, there are astrophysical considerations too that constrain axion 
coupling to photons or fermions. For instance, if they exist, being very weakly interacting particle, 
axions can drain away energy from stellar interiors. Since they are produced through processes like, 
e+ -b e” —i- 7 -b a or the cross channel reaction, e” -b 7 —> e” -b a or ^piasmon —>7 + a , j j ^ a 
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etc.; for a given coupling constant and mass, one can estimate the rate at which the axions would draw 
away energy form the stars. The bounds are placed by fact that, for a given energy budget of a star, 
the amount of energy drained by axion emission should be less than it’s observed luminosity. Apart 
from the ones mentioned above, there also have been experimental search for solar axions, through 
the conversion of an axion into a photon in a cavity, in presence of an external magnetic field. These 
searches also have placed some bound on the axion photon coupling. Incidentally its worth noting 
that since most of the astrophysical objects are associated with magnetic field the same process (or the 
reverse of it) can take place even in astrophysical environments too. 

The coupling of axion to photon is realized through a term in the Lagrangian of the following form, 

C=4-a^-B. ( 1 . 2 ) 

Where a is the axion and M is the axion coupling mass scale. The experimental bound on M coming 
from the study of solar axions is set to be, M > 1.7 x [8]. It may be worth pointing out here 

that, though it is usually believed that, M > 1.7 x lO^^^^Gel^, but this is a model dependent number. 
It should be noted here that the PQ symmetry breaking scale parameter is proportional to M. A 
detailed survey of various astrophysical bounds on the parameters of axion models and constraints on 
them, can be found in Ref. [9]. 

Since most of the bounds on axion parameters arise from astrophysical and cosmological studies where 
medium and a magnetic field are present, it becomes important to seek the modification of the axion cou¬ 
pling to photon, in presence of a medium or magnetic field or both. Particularly in some astrophysical 
situations where the magnetic component, along with medium (usually referred as magnetized medium) 
dominates. Examples being, the Active Galactic Nuclei (AGN), Quasars, Supernova, the Coalescing 
Neutron Stars or Nascent Neutron Stars, etc., to name a few. Apart from the ones discussed before, 
lately there seem to be some observational signature for possible existence of astrophysical objects,( 
called Magnetars ), with magnetic held strength, B ~ 10^® — 10^^ G, i.e., signihcantly above the critical, 
Schwinger value Be = rnlle ~ 4.41 x 10^^ G [10,11]. However we would like to emphasize that even 
normal astrophysical objects are always associated with magnetic held, though the strength of the held 
may not be as strong as Be- Thus justifying the role of magnetized medium on axion properties. Also 
recently there has been an attempt to describe the observed faintness of Type la Supernova remnants, 
based on axion photon oscillation in the magnetic held of the intergalactic medium. 

In view of these interesting physical applications of axion physics in astrophysics as well as cosmology it 
seems timely to hnd out the effect of medium and magnetic held, on the couplings of axions to photons. 

In this note we would investigate the matter induced photon axion coupling in a magnetized medium. 
Where the particle in the plasma will be considered to be mostly electrons, though it could be any 
fermion. Since the temperature in these astrophysical objects are not too large (of the order of hundred 
MeV or so at the most) this seems to be a reasonable approximation. Of course our formulas are general 
enough to be extended to any temperature and density. 

The organization of this document is as follows, in section II we would discuss about the physics of 
axion photon coupling and the model dependent uncertainties that enter in the axion photon coupling 
parameter M given in eqn. [1.2]. In the same section we would also try to give a brief over view of the 
existing axion models and their type of coupling with fermions. In the next section (i.e. section II), 
following Schwinger’s approach [12], we would elaborate on the details of the magnetized propagators. 
As would be discussed later, matter induced Axion photon coupling in a magnetized medium has two 
contributions in it, one coming from the magnetized vacuum and the other from the magnetized medium. 
Sections III and IV would deal with the details of those contributions. Finally at the end we would 
conclude by justifying our results through a general analysis and order of magnitude estimation of the 
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modifications one is getting from the presence of magnetized plasma. Finally we would like to conclude 
by pointing out the possible applications of our result. 

II. AXION MODELS. 

After the brief motivation of section [I], we would like to review the relevant details for the popular ax- 
ion models, since they would be useful for the medium induced modifications to axion photon coupling. 
Usually there are three types of invisible axion models found in the literature, (a)the Dine-Fischler- 
Srednicki-Zhitnitskii, usually referred as DFSZ [15] model, (b) The Kim-Shifman-Vainshtein-Zakharov 
(KSVZ) model [16] and lastly (c) the variant invisible axion model (VIA) [17]. The DFSZ axion has two 
doublets, (j)i{i = 1, 2) and one singlet y Higgs fields; the KSVZ model contains one (j) and one y Higgs 
fields along with a super heavy exotic quark that is singlet under SU(2)xU(l). The VIA model has 
some similarity to that of DFSZ model, except for the fact that it carries an extra Higgs singlet whose 
phase is identified with the axion. The most widely discussed axion model in the literature are DFSZ 
and KSVZ models. Apart form the what was mentioned at the beginning, the main difference between 
them comes out on the basis of their tree level couplings with leptons and quarks. The KSVZ model 
does not have any coupling to electrons at the tree level, however the same might be generated via 
radiative correction through photon-photon-axion vertex [18]; and is higher order in coupling constant 
hence not interesting. 

The U{l)pQ transformation rules for the Higgs fields DFSZ model are given by [15], 

(/)„ ^ (t)d and y ^ (H.3) 

Where u and d are the generation indices and couples only to up type quarks and similarly The 
Higgs field y interacts with and (j>d through the potential term. The transformation laws for the 
fermions under the same transformation are fixed by demanding that the Yukawa interactions would 
remain invariant. In the passing it may noted that, following [19], one can assume for convenience that 
the left handed doublets transform like singlets under U{l)pQ, thus fixing the PQ charges of the right 
handed fermions through their coupling to the specifically assigned Higgs field. 

Electromagnetic coupling of axions to photons were derived by Kaplan [20] and Sredeniki [21], [22], 
using current algebra techniques. In order to make the PQ current color anomaly free, a linear combina¬ 
tion of chiral currents for lights quark are usually subtracted from the same to define the new anomaly 
free axial vector current, given by, 

{uipl5U + zd'y^-f5d) . (H.4) 

Here, Apg is the color anomaly of the PQ charge, defined in terms of the generators of the group SU (3)c 
i.e., Xa (a=I to 8) in the following way, Sab^pQ = Tr[AaAbXf]. Where, the trace is understood to be 
taken over all the Weyl fermions and Xf are their generation specific PQ charge. 

We note for the sake of completeness that, four divergence of Eqn. [H.4] yields, 

= TT[XfQ}]-A‘pQ^^^^^ -^^(2tmuU^bU + 2ztmddjbd). (H.5) 

The Lagrangian describing the axion-fermion ( to be considered as lepton for the estimates made in 
this note ) interaction is given by [23], 

^af = ^ (HJ) 

Ja . 
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where g'af = ^f'T^f/fa is a dimensionless Yukawa coupling constant, XfS,{ the model-dependent factors 
) are the PQ charges of different generations of quarks and leptons [9], as given in eqn. [II.4]. Lastly 
TO/ is the fermion’s mass. In eqn. [II. 6 ] the sum over f stands for sum over all the fermions, from each 
family. Although, in places instead of using the Lagrangian given by, [II. 6 ], the following form for axion 
fermion Lagrangian has also been employed, 

/:a/ = -2*5l/^(§/75VI'/)a, (11.7) 

/ 

but, the correctness using [11.6] has been pointed out by Raffelt and Seckel in [24]. For our calculations 
we would however use [ 11 . 6 ], with gaf = 


The necessity for going through the details of the model becomes more apparent, if one recalls the way 
bounds were derived for M, in eqn. [1.2]. The use of the form of that Lagrangian leaves much scope for 
the uncertainty in estimating the actual PQ symmetry breaking scale — coming from model dependent 
assignments of PQ charges of the fermion multiplets. To elaborate this, point we note that, the effective 
photon photon axion coupling is obtained via the electromagnetic anomaly generated in the processes of 
making the PQ current free of anomalous divergence coming from the strong interaction sector. Thus, 
the axion photon photon Lagrangian as obtained from eqn. [11.5] [23], is given by: 


C 


a-yy 


327r2/a 


A em 

^PQ 


C 2(4+ z) 
^'^3(1 -kz) 


aFF. 


( 11 . 8 ) 


In eqn. [11.8], following are the definition employed, A|lg = Tr[Qf • Xf] and z is defined to be the ratio 
of the masses of two light quarks, i.e z = ^ and Qf is the fermion electric charge, in units of electronic 
charge. The axion photon mixing Lagrangian in an external magnetic field turns out to be. 


r = - 




327r2 


aFF 


Ext 


When ga-^^ = — 

la 


Aem _ AC 


2(4-kz) 


PQ ^pQ3(TT^J 


(11.9) 


The mixing strength, M as had already been given in eqn. [1.2] turns out to be, M = —. It is now 

rather easy to see that as one takes into account the effect of matter contribution (as shown in Fig.[l]), 
there will be additional contributions to 5077 given in [11.9] and hence to M. Though in principle (with 
out matter effects) the choice of PQ charges can make coupling constant 5077 extremely small, but as 
we would show in this paper that as one takes into account the magnetized matter effects, there are 
additional modifications to the vertex, generating momentum dependent coupling in addition to the 
details of the model and the details of the PQ charge assignments. We would try to establish the same 
in this note. 


III. AXION PHOTON VERTEX IN A MAGNETIZED MEDIUM. 

A. Fermion propagator in a magnetized medium. 

We are interested in physical processes in an external background magnetic field. Without any loss of 
generality, the same is taken to be in the z-direction and would be denoted as B. Any charged fermion 
propagator in such an external magnetic field, in Schwinger’s approach [12-14] is given by: 

POO 

iS^{p)= dse‘^‘'P’^'>G{p,s), (III. 10) 

Jo 

with $ and G defined in the following way: 
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(III.ll) 


«.(p,.) s is (p; - pi - m») - e|.|, 

^ieQ fBscTz / ^—ieQfBscTz 

" cos{eQfBs) V " cos{eQfBs) 

= \[l + iazta\i{eQfBs)){i)\\+m)+sec^{eQfBs)i>A , 


ll>± 


(III.12) 


The quantity Qf stands for the charge of the respective fermions, in units of electronic charge. Also 
the following useful relation should be noted, 


CTz = illl 2 = -707375 , 


(III.13) 


and we have used, 

^reQjBsass ^ cos{eQfBs) + i(T, sin(eQ/Ss). (111.14) 

It should be noted that, the metric convention for this propagator is To be more specific, 

according to the convention we follow, 

i>\\= "1013° + 73P° 

i>± =llP^ + 12P^ 

2 2 2 

P\\ =P0-P3 

2 2,2 

P_L — Pi + 1 ' 2 - 

Of course in the range of integration indicated in Eq. (III. 10) s is always positive and hence |s| equals 
s and by virtue of the e prescription the exponent damps out at s equal to infinity. Following standard 
prescriptions of thermal field theory, in the presence of a background medium, the above propagator is 
modified to [25]: 

iS{p) = iS^{p) + Sl{p), (III. 15) 

where 

Slip) = -Pf(p) 

and 

sl{p) = joSlHp)'to, (111.17) 

for a fermion propagator, such that 

/ OO 

dse^(P>^)G(p,s). (III.18) 

-00 

The information about the medium is carried by, ?7 f(p) that in turn carries the information about the 
distribution function for the fermions and the anti-fermions: 


iSlip) - iS^ip) 


(III.16) 


?7f(p) = 0(p-w)/f(p,P,/3) + &{-p- u)fFi-p,-p,P), (III.19) 

fp denotes the Fermi-Dirac distribution function: 

jFip.isJ) = ^ ^ . (IIL20) 

and 0 is the step function given by: 

<d{x) = 1, for X > 0 , 

= 0, for X < 0 . 

Here the four velocity of the medium is u, in the rest frame it is given by = (1,0, 0,0). 
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IV. EXPRESSION FOR PHOTON AXION VERTEX IN PRESENCE OF UNIFORM 
BACKGROUND MAGNETIC FIELD AND MATERIAL MEDIUM. 

In order to estimate the loop induced 7 — a coupling, we would start with the Lagrangian given by 
Eqn. [11.6]. The effective vertex for the 7 — 0 coupling can be written as ( for the sake of brevity, we 
define the notation, p' = p + k): 

^^{k) = {-igafeQf)kf^j ^^^Tr [7^75iS'(p)7^tS'(p')] • (IV.21) 

The effective vertex given by [IV.21], is computed from the diagram given in [Fig.l]. One can easily 
recognize that, eqn. [IV.21], has the following structure, T„{k) = k^Il^j,{k). Where is the axial 
polarization tensor, comes from the axial coupling of the axions to the leptons and it is: 

(-l)(-*)^(5a/eQ/)y^|^Tr[7^75iS'(p)7,,tS'(p')] • (IV.22) 

The axial polarization tensor, 11^^^, would in general have contributions from pure magnetic field back¬ 
ground, as well as magnetic field plus medium (in the text, the same might as well be referred as 
magnetized medium). The Pure magnetic field contribution (i.e the contribution devoid of any thermal 
phase space factors) and the one with magnetized medium effects, are to be found in the following 
expression. 


+ l^llb^S'^{p)-i^S]^{p')\ . (IV.23) 

The pure magnetic field contribution has already been estimated in [26-30]; however we have rechecked 
and verified the same according to our convention and would report about it in the next section. 
Following that, the thermal part would be the one, we would be dealing with. Incidentally since we are 
dealing with the dispersive effects, the contributions coming from the absorptive part of the diagram 
are ignored (i.e we have ignored the appearance of the terms in Tr ['^^i'ybS'^(p)')vS'^(p')] in side the loop 
integral ). This is justified further by the fact that we are interested in photons far below the threshold 
of pair production. 


A. Contribution From Magnetized Vacuum. 


The VA response function in a magnetic field has been calculated in Ref. [26-30]. However since 
our conventions are different; also since there are some discrepancies in the results stated there ( mostly 
typographical though ), we have redone the calculations and our result reads, 


KHk) = 


(47r)2 


^00 ^-1-1 

/ dt dv e 


4>o 


{( 


1 — 


k\\ 2TOg^ F]u 


(Fk) 


+ R[k^^{kF)^ + k^^ikF),'\'}, (IV.24) 


Where, R 


I — n sin Zv sin Z — cos Z cos Zv 
sin^ Z 


and (j)o =it 



cos vZ — cos Z 
2Z sinZ 



(IV.25) 


and with ^ is the dual of the field-strength tensor, with Z = eQ fBt. Therefore, 

following Fqn. [IV.21], The photon axion vertex in a purely magnetized vacuum, would be given by, 
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(IV.26) 


T''{k) = k>^U^-{k) 

The expression for r'^(fc) turns out to be, 


T''{k) = 


WafieQff 


(47r)2 

+ R 


pOO p-\-l 

/ dt 

lo J - 


dv e‘^° I ^ 


1 — 


kf\ — 2mi 




(kF) 


^ {kF)i, , 


(IV.27) 


This result is not gauge invariant. However, one may integrate the first term under the integral by 
parts [28,29] 


/ — 7 ^ —fcij — 2 to^ j = 2i+ dtk‘]_e'^°R. 

Jo ^ 2, J Jq 


(IV.28) 


to arrive at. 


r"(fc) = 


igafjeQff 


4i{kF)^ + dtj^ dve’^'>{(^-{l-v^)kfj{kF)^ + R k^^k^^{kF)^ 


(47r)2 

_igaf{eQff 


(47r) 


l‘00 /. + ! 

4i — / dt 


/o J- 


dve’^'>(l-v^'jkl 


(kF), 


(IV.29) 


As is evident, terms proportional to R in eqn. (IV.29) cancels out against each other leaving the vertex 
gauge invariant. Incidentally, in our way of writing (kFY = k^^F^'^ and {FkY = F'^^kf^. So the 
Effective Lagrangian giving rise to axion photon coupling in a magnetized medinm would therefore be 
given by, 


= aA"r,(A:) (IV.30) 

It should be noted that, a in Eqn.[IV.30], denotes the axion field. 

In the limit of w << m/, axion photon vertex can be written as: 

fci 4 ' 

4 H —— — 
eB 6h\ 

Y \ ■ 

^ J aF^^F^r (IV.31) 

Equation [IV.31], describes the effective axion photon interaction. 


-1 




-1 

327r2 


gafYQf) 


“ + 1 : 


V. CONTRIBUTION FROM THE MAGNETIZED MEDIUM 

Having estimated the effective axion photon vertex in a purely magnetic environment, we would now 
estimate the same in a magnetized medium. As before, one can do that by using the form of the fermion 
propagator in a magnetic field in presence of a thermal medium, as given by expressions(HI.IO) and 
(HI.18); on doing that, we arrive at: 
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FIG. 1. One-loop diagram for the effective axion electromagnetic vertex . 


/ j4 /'Cx> 

pOO 

/ ds'e'^^P ^ [Tr [ 7 ;, 75 G(p, s)j„G(j)', s')] Vrip) 
Jo 

+ Tr [l),75Gi-p', s')j^G{-p, s)] t]f{-p)] 
= -i-ifi9afeQf){-l) J J dse^^P’^'l 

pOO 

X / ds'e^*^^'’'*')R^^(p,p',s,s') 

Jo 


(V.32) 


where R^iy(p,p', s, s') contains the trace part. It should be noted that, in general R^jy(p,p', s, s') is a 
polynomial in powers of the external magnetic field ( with even and odd powers of B), i.e, 


RAi>^(p,p',s,s') = {p,p', s, s') + TZ^^\p,p', s, s') 


(V.33) 


In the sections below we would evaluate them one after another. 


A. R^ii/ to even and odd orders in magnetic field 

We denote the pieces even and odd orders in the external magnetic field in R^,y as R^^^ and R^^^ 
In addition to being just even and odd in powers of eQfB, they also are odd and even in powers 
of the chemical potential (the same would be clear as we go along) i.e., under charge conjugation, 
BSzp (—/r)&(—6), they both behave differently. More over their parity structures are also different. 

These properties may come very useful while analyzing, the structure of axion photon coupling, using 
discrete symmetry arguments; and that is the reason, why we have treated them separately. Calculating 
the traces we obtain (details for is provided in the appendix, in any case on contraction with 
vanish and if that can be established on general grounds, one does not have to evaluate it.), 

=4f?7_(p) [e^i,„||/3||p“iip''^ii(I-Htan(eQ/;Bs)tan(eQ/Ss')) 

+ sec^(eSs') sec^{eQfBs) 

+ sec^{eQfBs) sec^(eQ/^s')] (V.34) 


and 


= 4*7+(p) [-m^e^i.i2(tan(e(5/^s) -f tan(eQ/Ss')) 

+ {9iia\\P°‘'''p'y^ -Fpi,a||P“"p(,^)sec^(eQ/Ss')|tan(eQ/Ss) 

+ -h5,.a||P'“"P;ii)sec^(eQ/;Bs)|tan(e(5/Ss') 
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(V.35) 






We would further express Eqn.[V.35] 

as a sum of two pieces, = 4 i? 7 _|_(p) [— 2 TO^£^^i 2 (tan(eQ/, 8 s) + tan(e( 5 /; 8 s'))] + R°^(p,p', s, s'). 

Where R°j,(p,p', s, s') is defined as follows, 

R^fj = 4 zr 7 +(p) [m^£:^i.i 2 (tan(e( 5 /, 8 s) + tan(eQ/, 8 s')) 




+ + 9ua^^P°‘"p'^Jsec^ieQfBs')'^tan{eQfBs) 


+ [i9^ia^^P'°‘"Pl'^^ - gti.yPa^^p'"^'''' + gya^^p'°‘"P^f^) 


+ {9^La^^p'°'"Pu^ -k5i.Q,iiP'“iip^^)sec^(eQ/,8s)|tan(e(5/,8s') 

(V.36) 

P+ip) = Pf(p) + Pf{-p) 

(V.37) 

V-{p) = Vf{p) - Vf{-p) 

(V.38) 


they carry the informations about the thermal nature of the medium. More over, it should be noted 
that,according to the convention followed in the text, 

5 ''ii = ao6^ + 036°. 

In the steps below it will be shown that, as Eqn.[V.36] is contracted with k^, it would vanish. Hence, 
only the non-vanishing contribution would be proportional to the integral over the four momentum and 
the parameters p, s and s' respectively of. 


^ — 2 m^ef_ii,i 2 {tan{eQfBs) + tan(eQ/, 8 s')) 


In the rest frame of the medium, p ■ u = pq. Thus the distribution function does not depend on the 
spatial components of p. In this case we can write the expressions of R^®j and R^°2 using the relations 
derived earlier [31] inside the integral sign, as 


dx E_ tan(eQ/gs') 

tan(eQ/,Bs) -f tan(eQ/,Bs') 

^ tanjeQfBs) 

tan{eQfBs) + tan{eQ jBs') 

taii{eQfBs)+tan{eQfBs')^ 

_ tan^(eQ/gs') 2 

tan{eQ fBs) + tan{eQfBs') 

tan{eQfBs)+taTi{eQfBs')^ 

_ tan^(eQ/gs) 

tan{eQ fBs) + tan{eQ fBs) 

= (*^ + by - sec^{eQfBs)p']_)^ 


(V.39) 

(V.40) 

(V.41) 

(V.42) 

(V.43) 


Coming back to eqn.[V.34] can further be simplified using Eqn.[V.39,V.40] and one arrives at. 
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=4ir7-(Po) [e^,.a||/ 3 ||P“"p'^"(l + tan(eQ/Ss)tan(e(5/;Bs')) + sec^{eQfBs') 

+ sec^{eQfBs)] (V.44) 

The last term in Eqn.[V.34]has vanished because, on using the identities as mentioned above,(i.e 
Eqn.[V.39,V.40], it turns out to be proportional to , hence zero. The = symbol sig¬ 

nifies that the relations above are not proper equations, the equality holds only inside the momentum 
integrals in Eq.(V.32). 


VI. EVALUATIONS IN EVEN AND ODD POWERS IN MAGNETIC FIELD. 
A. Evaluation of to even orders in the external field 


The axial polarization tensor, even in powers of the external field ( denoting by is given 


by: 


n 




= -{-^en- 


pOO ^OO 

J ds' (VI.45) 


Using Eq.(V.44) in the rest frame of the medium, we have: 

= 4ir7_(po) [£^J.ua^^ 0 ^^p“"p''^" (1 + tan{eQfBs) ta,n{eQfBs')) 

+ sec^{eQfBs') + e^va^ 0 ^^p°‘^p^'''' sec^{eQfBs)] . 


(VI.46) 


Noting that it is possible to write, q°‘Pa = q°'''Pa« + q°‘^Pa±_ eq.(VI.46) can be written as. 


=4ir7-(4'o) [{£tiuapp°'p'^ - £,j.va 0 ^P°‘p'^^ - tan(e(5/.Bs) tan(eQ/.8s')) 

+ £p.uap^p“p'^^ sec^(eQfBs') + £^„a^ 0 P°^p'^ sec^ieQfBs)] . (VI.47) 

Here throughout we have omitted terms such as £ij.i^a_i_P±P°‘^P'^^ ^ since by the application of Eq.(V.39) 
we have 






ta.-n{eQ fBs') 

tan{eQfBs') + tan{eQfBs') 


which is zero. 

After rearranging the terms appearing in Eq.(VI.47), and by the application of Eqs.(V.39) and (V.40) 
we arrive at the expression 


= 4ir?-(po) 


^fiiya^P k^{l + tan{eQfBs)tan{eQfBs')) + £fj,uapj_k°‘k^-^ 
, taa{eQfBs) — tan{eQfBs') 


X taii{eQfBs) taii{eQfBs' 


tan^eQfBs) + taTi{eQ fBs') 


(VI.48) 


Because of the presence of terms like e^uapk^ and £^uap_i_k°‘ as we contract by k", it vanishes. In 
view of this result it is tempting to conclude that, axion coupling to photon, in a magnetized medium, 
cannot be even in powers of external magnetic field. With all possibility this would be forbidden by 
the discrete symmetries i.e CPT. However, apart from the discrete symmetry arguments, there is an 
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alternative, recalling the fact that each power of eQfB actually denotes insertion of very soft photon 
(i.e all components of the four vector, ^ 0 ). In this case, the coupling of the dynamic photon along 
with the axion to even number of soft photon insertions, makes the power of the electromagnetic vertex 
odd. Since axion has spin zero, in order to match the total spin of the system, the sum over photon 
spins should add up to zero, which is impossible with odd number of photons. Hence this term should 
vanish in principle. Which is the result one arrives at after explicit computation. 


B. The Result for Odd Orders In Field Strength 


As has already been noted that, to odd orders in field strength one has and it can be expressed 
as a sum over two terms, one proportional to m? and the other is given by. 


=4ir7+(p) m^£^^i 2 (tan(eQ/Hs) +tan(e(5/Hs')) + | 


+ +g^Q,iiP“iip'^^)sec^(e(5/Hs')}tan(eQ/Hs) 

+ - gtiuPa^^p"^" +5i^a||P'“"PM||) + + 5i^a|| ) SeC^ (eQ/^s) | 

X tan{eQfBs')]. (VI.49) 


In this subsection, we would outline the proof that Eqn. [VI.49] on contraction with vanishes, for all 
v. The proof follows in two steps, first one demonstrates that for = 1 or 2 it goes to zero. Followed 
by that one shows that the same vanishes for i/ = 0, or 3 as well. Denoting ^ = 1 or = 2 collectively 
a,s , v±, we show in appendix B, 

k^Ri°}^{k) = 0. (VI.50) 

Next we claim that, for ^ = 0 or 3, denoted collectively as (i.e the longitudinal components); 
would vanish. One can verify, with a bit of algebra, that the same can be written as a sum over two 
pieces. That is. 


d(o) _ jLM/3 (o 1) r-Ai E?(o.2) 

Apart from the uninteresting overall constants and integrations over the variables p,p’, s and s’ the 
pieces are found to be proportional to. 


- P;p + 2(A: -p)^ (tan(eQ/Hs) + ta,n{eQfBs')) 


k^{ta,n{eQfBs) — tan{eQ fBs')) 


(VI.51) 


The fact that this piece vanishes exactly on general symmetry grounds of the integrals was first noted in 
[33] and also in [32], so we skip the details till appendix C. The required elaboration is provided there. 

The left over pieces, not considered yet, is That is given by. 


kf'R 


( 0 . 2 ) 

fiO 


—4f?7_|_(p)fc^ 



(tan(eQ/Hs) +tan{eQfBs')) — k^tan^{eQfBs') 
sec^ eQfBs 

tan{eQ fBs) + tan{eQfBs') 


Similarly, 


(VI.52) 
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= 4zr7+(p)fco 



{taii{eQ fBs) + tan{eQfBs'j) — k‘]_taia^{eQ jBs') 
sec^ eQfBs 

tan{eQfBs) + t&n{eQfBs') 


(VL53) 


In appendix D we have shown that Eqns [VL52] and [VI.53] vanishes.Thus the proof that, 
, s, s') = 0 is complete. 


VII. TO THE EFFECTIVE INTERACTION. 


So using the results of the previous section, the effective axion photon vertex can be written down. 
Noting that, 

= -Siw?r]+{p) [fc'"£:^i,i2(tan(e(5/.8s) + tan(eQ/,Bs'))]) (VII.54) 

the vertex function Ty{k) ( using Eqn. [V.32]) turns out to be, 

^uik) = {gafsQf) J-^^r]+{p) j ds Jds' 

X [tan(eQ/^s) + tan(e(5/.8s')] (VII.55) 

Since the perpendicular components of the momentum integrals are Gaussian, one can carry them out 
without any difficulty with the following result, 


r^{k) = (gafeQf) (8m^fc^e^^i2 


(Bp 

(27r) 


^OO pOO 


g+{p) / ds / ds'e 

J — oo J 0 




ta,n(eQ j: Bs) ta,n(eQ j'Bs') 

^ -ieQfB[tan{eQfBs)+tan{eQfBs')] ^-^eQfB [t,„(eQ,B3)+tan(e0f bX)] 

47r [tan{eQfBs) + tan(eQ/.Bs')] 

(VII.56) 


Exact evaluation of [VII.56] is extremely difficult, however it is possible to get some analytical results 
under certain approximations. In the long wavelength limit and for oj < mf one can get an analytical 
estimate of the same. However for energies above the Pair Production Threshold (PPT), the one has 
to be very careful in evaluating the vertex functions. In the passing we note that, in the limit of 
m > p Eqn. [VII.56], can be evaluated analytically using Poisson summation formula and Laguerre 
polynomials. However in the limit of vanishing |A:_l|, the leading order result can easily be recovered 
from eqn. [VH.56] by setting setting |fcj_| ^ 0 in the same. The issue of evaluation of the same for all 
values of to and k would be taken up in a separate publication. Since, in this work, we would like to 
focus our attention to photon energies below pair production threshold. In what follows, from now on, 
we assume photon energy to be below PPT and would take the long wave length limit, i.e [fej.] ^ 0, 
and try to estimate the magnitude of the contribution of the vertex function. 


lim r^(fc) 
IfciHo 


(gafeQf) (8m^A:^£^^i2) 


-ieQfB /'d^Pii 


^oo poo 


Att j ( 27r )^ 
d?p 


77+(p) ds ds'e 

J — ooJ 0 


/ is{p'^ ——m^)) 


-(igaf(eQf)^) (4m2fc'"e^^i2) B J j^p+(p)S{p‘^^ - w?)) 


(p/2 _ 


(VH.57) 


The last step in Eqn.[VH.57]follows from the definition of the delta function preceded by subsequent 
integration over s'. 
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(VII.58) 


/ oo 

-oo 


is{p'^ —m^) 


ds = 27r5(p^ 


r*oo 


e II ds = 


—i 


bf - 


— m^)) 


m?) 


In the rest frame of the medium p.u = po, in this frame, the thermal factor, pp{p.u) turns out to be, 
p+ipo) = npilpol^d-) + npilpol, —fJ-) [33], (Here npix), the thermal distribution function, is defined as 
npixjp) = , as is the case for Fermi distribution fn.). using the same and the delta function 

constraint, one can simplify Eqn. [VIL57] further, to arrive at. 


r^{k) =-{gaf{eQff) 


d'^p,, 


nF{\po\,p) +nF{\po\,-p) 


(27r)2 L 

y.5{pl- m^) 


1 


(P + 2(p-fc) ) 


\2\ / 


= -16(5a/(eg/) ) Hk^ ,k-u,P,p) 


(VII.59) 


Where A(fcj^, k^^ ■ u, (3, p) is, 

A(fci^,fc • u,P,p)=J<fp^^ 


nF{\po\,p) +nF{\po\,-p) 


rrrS{p^ — m^) 


(fc2 + 2(p.fc)„)y 

Therefore the effective Lagrangian for axion photon interaction would be given by. 

F F^^ 

jrBF,P = = -16 (gafieQff) ^|^A(y, k-u,(3,p). 


(VII.60) 


(VII.61) 


Equation [VIL61] is the expression for axion photon coupling in a magnetized medium, in the limit, 
jfcxj ^ 0. It is important to note that, this equation depends on many physical parameters, e.g., the 
temperature of the medium ( /? = 3/T), number density of the fermions (which in turn is related to p), 
mass of the particles in the loop (m), energy and longitudinal momentum of the photon ( i.e. fc||) and 
of course the symmetry breaking scale ( which is included in the coupling constant, gaf)- 


Since the basic purpose of this exercise is to find out the effect of the medium on the mixing of axions 
to photons, we would have to choose some energy scale for the parameters entering into the statistical 
factors. Since most of the studies in literature deals with optical photons, we would assume the following 
scale, i.e TO/ ^ and finally to > ks. We would like to point out that, to = ko and ks would henceforth 
be denoted as k. Armed with these we can now turn to evaluate, Eqn. [VII.60] .In the limit as mentioned 
above. 


A{k‘^^,k-u,f3,p) 


POO 

/ dp 

nF{Ep,p) + nF{,Ep,-p) 

2Ep 

/ 

/ —OO 





(VII.62) 


We would like to point out that, while coming to the last step in equation [VII.62], we first had 

p 2 p.k—k^^ 1 

performed the the po, integral, then had expanded the resulting denominator, in powers of 




while retaining the leading order term in powers of 


r 


UjE-n 


. In the expressions above Ep = a/( p^ + to^). 


Hk^ k-u,P,p) = - 


dp 


nF{Ep,p) + nF{,Ep,-p) 


2Er, 


\2Ef 


(VII.63) 
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nF{Ep,fj.)+nF(Ep, — /i) 
2 Ep 


From now on we denote A = — J^^dp 

its easy to see that A(fc^, k ■ u, (3, p) = A = — dp 
small A:&fc||), A(fcj^, k ■ u,l3,p) = (77) 


If we neglect pieces 

nF(Ep,iJ.)+nF(Ep,-fi)'\ rn^ 

2Ep 2El ■ 


proportional to 
Other wise (for 


A. Limiting Case: m S> p. 

Analytical evaluation of eqn. [VIL63] for arbitrary values of chemical potential p and temperature 
is an extremely difhcult task. However as has already been mentioned, under some approximations 
it may be possible to get an analytical form for the same. One such approximation is m ^ p. This 
is reasonable in most of the astrophysical and all of cosmological situations, barring the core of the 
Supernova or active galactic nuclei where this approximation may not hold. In this approximation,( i.e 
p) thermal Fermi-Dirac factors can be expanded in a series to give. 


npiEp, p) + np^Ep, -p) = 2^(-I)” cosh [{n + l)l3p] e 


(VII.64) 


n—0 


In order to evaluate, Eqn. [VII.63], we would substitute Eqn. [VII.64] in the same eqn., to arrive at, 

^-{n+l)l3Ep ^2 

—1)'‘cosh [(n + l)/3^] j d} 

n—0 


00 « 

A = — ^(—I)" cosh [(n + l)/3^] dp- 

^—n 


E, 


p 


E'^ ' 

p 


(VII.65) 


The evaluation of Eqn.[VII.65] is a bit difficult, however, expression appears a bit less formidable, if 
one identifies, ^/s with Ep and uses the following integral transform. 


S 2 e “'Z® = TT 2 


si du 


(VII.66) 


to convert Eqn. [VII.65] close to a Gaussian. Upon using this transform, the p dependent part of the 
integrand goes over to, 


dp 


o-p “ 


p^ + 


= — [l — ^{m\/u)\ e“ 


m 


(VII.67) 


Here the function ^{my/u) stands for error function. Fortunately there exists another integral repre¬ 
sentation for the right hand side in Eqn. [VH.67], given by. 


m J m 


(VH.68) 


Use of this expression makes the expressions manageable. Using eqn.[VH.68], and expanding the left 
hand side in powers of u and then performing the t and u integration one arrives at. 


A =+DM f: 




n—0 


1^0 


i\ 




3"* t'^dt. 


(VH.69) 


dllt should noted that this s is different from the proper time parameter s introduced in the fermionic Schwinger 
propagators. 
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In Eqn.[VIL69], u integration would result in giving modified Bessel function and following that, the t 
integration, can be performed using the formula provided below, 


P 


(;)■ 


to arrive at the following expression for A. 


A = -I cosh|(„ + 1)M f; 


n—0 


Z=0 


i\ 


, (n + l)/3m, i+i 

{- -= K{i+i) \(n + l)m0\. 


(VII.70) 


(VII.71) 


One can use the expression [VII.71] in equation [VIL61] and estimate the contribution of magnetized 
medium to photon axion mixing. For simple cases, e.g. in high temperature (more relevant for the cos¬ 
mological and astrophysical scenario) or low temperature (intergalactic regions) it should be possible 
to use different expansions for the modified Bessel functions and estimate the size of the corrections to 
arrive at the analytical results. 


B. Limit: m < fi < oo 


In certain physical situations matter density can become extremely high, i.e the highly degenerate 
Fermi system. Examples are nascent neutron stars or the core of the core collapse supernova where 
matter density can be 10^^^ or more and temperature can vary between several tens of MeV to few 
KeV. Keeping this physical situation in mind we would try to estimate the contribution of A for various 
values of chemical potential. Noting that the 07 effective Lagrangian can be written in the form 




oFF 

647r^ 


9af(s.Qf) 



oFF 

327r2 


9af{sQf) 



16-A 


1 


g/3m(p2 + i)^_/3^ _|_ 2^ 


+ - 9 ) 


(VII.72) 


The integral inside square bracket in Eqn. [VII.72 ] is being denoted by A where A = —2A. It is 
worth noting that inside the integral, the function inside the braces for m <C /i and limitT ^ 0 can be 
approximated by 0(^ — p). In this limit we can write, 


Ltx—►0-^ 


9- 

dp 

Jo (p2 -f 1)3/2 


= sin 


tan 3 




(VII.73) 


In the limit jj, m, the right hand side of Eqn. (VII.73) ~ 1. 

Eqn. [VII.72] has also been estimated numerically, for m < /i < 6 when the temperature is held fixed 
at, T = 10“3 in units of fermion mass. The behavior of the function A in this range can be seen in 
[Fig.2 ] which saturates at 0.98 providing good agreement between analytical and numerical result. 


The important part of this analysis is, for sufficiently large values of matter density and low temperature 
the value of A is seen to saturate at a value not exceeding one. Though for smaller values of matter 
density it ranges between half to unity. 
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VIII. DISCUSSION AND CONCLUSION 


As has already been mentioned at the beginning, the purpose of this note is to find out the modifica¬ 
tions in the axion photon vertex in the presence of an external magnetic field and matter. Field induced 
part of the photon axion vertex has been worked out earlier in [34], however the effect of magnetized 
matter wasn’t taken into account. In this note we have considered the effect of both external magnetic 
field and medium on the vertex. More over we have also written down the explicit formula for the 
vertex involving the PQ charges of the fermions. In the light of these estimates, it is possible to write 
down the axion photon mixing Lagrangian, for low frequency photons in an external magnetic field, in 
the following way: 


(viii.74) 

Where is the usual axion photon mixing term in vacuum, resulting form the electromagnetic 

anomaly, and is given by [23], 


L 


vac 

a-y 


c 


B 

ay 


rB,ti,0 

^ya 



2{4 + z) 
3(1 -f z) 


(VIII.75) 


To remind ourselves, in the equations above, z = where niu and stands for the masses of 
the light quarks. The anomaly factors are given by the following relations, = Tr{Q‘j)Xf and 
= TrXaXbXf (where the trace is over the fermion species.). Hence, in the limit of [fcj,] —> 0 and 
to << ruf, using eqn. [VIII.75], one can write the total axion photon effective Lagrangian: 


L 


T otal 
ay 


gaj'i + 4 -F 



-16 


f 


327r2 


oFF. (VIII.76) 


In order to describe the axion photon interaction, for static magnetic field and real photon, the effective 
Lagrangian employed in the literature is usually given by the following relation, 

—aE.H“‘. (VIIL77) 

Where the parameter M, defines an energy scale in terms of inverse of the axion photon photon coupling 

oo 2 

constant in vacuum, i.e M = —. In principle this factor is a model dependent quantity, and 

depending on the particular choice of PQ charges, 5077 can vary between zero to hundred. This is one 
of the model dependent uncertainty. On the other hand, in some experiments, attempts have been made 
to give a bound on M thus ga-y-y through the detection of solar axions in a magnetic cavity haloscope. 
The astrophysical bounds on M come from the constraint that the stars don’t lose energy through axion 
photon conversion at a rate faster than the generation of energy. As a result of these investigations, 
the present day bound on the energy scale is believed to be, M < lO^^GeU. On the other hand in 
this note, we have been able to point out that apart from the model dependent uncertainties, there 
is also medium dependent uncertainties that can introduce some variation on the bounds on M. The 
amount of this variation would depend on the kind of environment one is interested in. The details of 
its implications for various models and physical situations are extremely interesting by their own right 
[35-44]; however these points are beyond the scope of this work. 
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T 


FIG. 2. A vs. for relativistic degenerate medium. The parameters are as follows: m < /i < 6 when the 
temperatnre is held hxed at, T = 10“® in units of fermion mass. 
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APPENDIX 


Appendix A : The Equations. 


Contribution from magnetized material part toVA response function comes through the pure thermal 
propagators (in the external magnetic field). Eqn.[IV.23] reads: 

(k) = j ^^^Tr [7^75S']^(p)7^fS'^(p')+7M75i5'^(p)7i^'S'B(p')] (H.78) 

Eqn.[H.78] has two pieces in it. In what follows we would try to demonstrate that except the thermal 
weight factors (i.e the ?7 f(p)s) the traces of both the terms come out to be the same. 


To begin with, we would concentrate on the first piece. In terms of the functions G(p, s) and G{p', s') 
the same can be written as. 


POO POO 

Tr [7^755*^(p)7i,iS'B(p')] =-Tr 7^75 / ds G{p, s)-f„ / ds' ’"'^Gip , s') 

J —00 0 


dF{p) 


pOO 

/ (is'e^(^’'’'*')Tr[7^75G(p, s)7,.G(p',s')]?7f(p) (H.79) 

> Jo 
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We would like to note the appearance of the integration variables,p, p' and s, s'. 


The second term in the right hand side of eqn.(H.78) needs some formal mathematical manipulations. 
Its better to do the same in two steps. The first step would be shift the integration variable p ^ p — k. 
And following this the next step would be to change the integration variable p —p. As a result of 
this manipulation one arrives at, 


Tr [- 1 ^,lbiS'^{p)luS''g{p')\ = Tr [-i^-i^iS'^{-p')-i^Sl{-p)\ 


= -Tr 


pOO pOO 

7/x75 / ds'e‘^^~P'’''''^G{-p',s')j^ / ds s) rjpi-p) 

1 / 0 J — 00 

/ OO POO 

/ g;g'g'&(p'.s')Xr[7^75G(-p',s')7i-G(-p, s)]pf(-p) (H.80) 

-00 J 0 


The complete expression is given by, 

*n;!r(fc) = -He2)(-l) J dse'^^P’^^ ds'e'^^P’’^'^ 

X [Tr [7 m75G(p, s)-f^G{p', s')] Pf(p) + Tr [7 m75G(-p', s')'y^G{-p, s)] pf(-p)] (H.81) 

=-{-ie^){-l) J J dse'^^P’"^ J ds' R^j,(p,p', s, s') (H.82) 


where R^i/(p,p', s, s') contains the trace part. It is worth noting that the two pieces appearing in 
the second line of Eqn. [H.82] are related to each other by the following transformation (p ^ —p') and 
(s ^ s'). So once we compute the trace of one of them, the other can easily be obtained by interchanging 
the variables,p ^ —p' and s ^ s'. 

Hence, in the following we would evaluate just one of them. If we introduce the following short hand 
notations. 


tan{eQfBs) = Ts and 

sec^{eQfBs) = Sg 


tan{eQ fBs') = Tgi and 

sec^{eQfBs') = Sg' 


TpS 


(H.83) 

The expression for, G{p, s) can further be written as. 

G{p,s) = {1 + iaztan{eQ fBs)){f>\\ + 

to) — sec^{eQfBs)f)± 

(H.84) 


= (I + ia;,Ts){f)i\ + to) - Ssf>j_ 


After doing a bit of algebra, the Dirac trace part of the first term of Eqn. [H.8I], odd in powers of 
eQfB turns out to be (written in terms of short-hand notations introduced in Eqn. [H.83]), 




d^p r 

( 2^)4 J-c 




X i 


{crzi>\\li>i)\^Ts +;i||7,.cr2/i'Tj^ -I- to^ {Tf^^az^iuTs + T^syi^a^Ts) 




(H.85) 


(H.86) 


In order to simplify things a bit, we note that, using the defn of az, (i.e CTz = — 7 o 7375 )one can write, 

= *7172/^11 = 7o73/^|| 75 = 7o73 (7oP° + IsP^) 75 = - (7oP^ + 73P°) 75 

= —7a||P^", (Q!|| could be, either 0 or 3). (H.87) 
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We would also like to point out that, depending on the value of an, would take its complimentary 
value, i.e to say, if an = 0 then aj^ = 3 etc. We also note in the passing that in this notation. 


= - 70,11 75P “ii . 

= - 7 o|| 75 p“" and finally, (H. 88 ) 

Tr[rp5cr2 7 i.] = -4e^i.i2 (H.89) 

Tr[r^ 57 ^CT 2 ] = -4e^i.i2. (H.90) 

Given these equations, it would be easier to evaluate the Dirac traces appearing in Eqn. [H. 86 ]. We 
would start with the piece proportional to rn^. It reads, 

= -4m^e^i.i2 + T'^). (H.91) 


It is worth noting that, the sign of the piece proportional to is negative, more over the piece is 
antisymmetric with respect to the indices /r and v. Next we would like to evaluate the trace of the 
pieces having parallel component of momentums p and p', i.e.: 


Tr 




= Tr 


7/x7q 



= 4 


+ 


(5/XO||P'“"PxX|| 


-5/xixPa||P'“" +5i^a||P'“"PAX||) T'^ 


(H.92) 


It should be noted that unlike the momentum dependent piece, this part is symmetric in the indices 
/i and V. In other words these terms are insensitive to the position of 75 . That is to say,inside the 
trace, if we would have pushed the 75 matrix to the p vertex, these pieces would remain insensitive. 
Its also important to notice that these pieces are symmetric with respect to exchange of the variables 
(p, s) (p', s'). The last trace that needs to be evaluated is. 


Tr 




5/xa„p“llp',x^ TsS: 


= 4 


+ ( gfJ.a,,p“"Pu^ +9 t'Ct, P°"PiJ.± Ts'Ss 


(H.93) 


As can be verified, Eqn. [H.93] is symmetric is the indices plzv, more over as before it remains the 
same under the transformation of variables (p, s) (p', s'). On the other hand the symmetry under 
the exchange of (p, s) {p',s') can be exploited to evaluate the second trace in Eqn.(H.81). Since 
the momentum dependent terms all appear in even powers of p and p',these equations are insensitive 
to the following transformation (p,p') (—P^ ~p)- Hence under the transformation 

{p, s} {-p', s'} , (H.94) 

the traces would remain invariant. Therefore the Dirac trace of both these pieces are the same and 
hence upon using the definitions of, Tg, Tgi, Sg and Sgi as given in Eqn. [H.83], we recover equation 
[V.35]. 


APPENDIX 

Appnedix B: Proof Of — 0 

The expression for R^°j^(A:) as follows from Eqn. [VI.49], turns out to be. 
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APPENDIX 

Appendix C: Proof Of = 0 

In the text, we claimed that the contribution coming from Eqn. [VI.51] must vanish. Here, we justify 
that claim. Vanishing of Eqn. [VI.51], actually amounts to showing that the following integral, 

/ 74„ POO POO 

V+ip) J + 2{k ■ p),,) (tan(eg/Ss) + tan(eg/6s')) 

— fc5_(tan(eg/,8s) — tan(eg/,8s')) =0. (H.99) 

We begin by defining two new parameters 

C= ^eQfB(s + s'), 

(:=^eQfB{s-s'). (H.lOO) 

and noting, ( following [31] ), that: 
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(yfc^ + 2 (p./c)|| - p'l sec^{^ - () + pI sec^{^ + O) ■ (H.lOl) 
Using Eqns. [V.41] and [V.42] in Eqn. [H.lOl] and little bit of algebra, one can further show that, 

/ t 4 „ pOO pOO J 

J ds'—J^{^,C), (H. 102 ) 

where 

= (tan eQ/, 8 s +tan (H.103) 

with s and s' related to ^ and ( through Eq. (H.lOO). We can now change the integration variables to 
^ and ( to arrive at (ignoring the unimportant constants in front), 

/ t 4 pOO pOO 1 

(p) j_^ d( 0 0 

/ t 4 pOO pOO r ri r 

= - / V+{P) £ , (H.104) 

since the other term vanishes at the limits. In this integrand, C = which means s' = 0. Looking back 
at the definition of iF, we find 

=exp|$(p, ^^^)|tan2^. (H.105) 

This is an even function of p, whereas ( p~ x ? 7 -|_ (p) I is odd. Thus, the expression vanishes on integration 

V "'ll J 

over p. 


APPENDIX 


Appendix D: Proof Of — 0 


In this appendix we show that, = 0. The two components of are 0 and 3. Respective 

components those need evaluation, are, 


= —4ip+(p)fc^ ^p^ — (tan(eQ/,8s) +ta.n{eQfBs')) — tan^{eQfBs') 


sec^ eQfBs 

tan{eQ fBs) + tan{eQfBs') 


(H.106) 


Similarly, 


= 4ir7+(p)fco ^p^ — {taTi{eQfBs) + taii{eQfBs')) — k^tan^(eQfBs') 

sec^ eQfBs 

tan{eQ fBs) + tan(eQfBs') 


(H.107) 


Expressions inside the square bracket of Eqns. [H.106] or [H.107] are the same, introducing compact 
notation. 
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Ii{p,k) = 4ir7_|_(p) {taii{eQfBs) + tan{eQfBs'j) — k^tan^{eQjBs') 

sec^ eQfBs 


tan{eQfBs) + tan(eQ/;Bs') J 


(H.108) 


Actually the expression defined by n(p, k), appears inside the integral over the loop momentum variable 
p as well as the parameters s and s'. Calling the integrals over p, s and s' of n(p, k) as n(p, k), the 
same is: 

d^p 


oOO pOO 


n(p, k) =1 


(27r)4 J_^Jo 

— k'^tan^ {eQfBs') 


dsds' Ui? 7 +(p) [(p^ - m ^^{tan{eQfBs) + tan{eQ fBs')) 
sec^ {eQfBs) 


tan{eQ fBs) + tan{eQ fBs') 


(H.109) 


Though not obvious from appearance, but we would show that Eqn. [H.109] vanishes identically. To 
show that, we note that it is possible to write. 


dse 


s(p -m - p^-e\s\^j2. _2 


(p^ — m^)(tan(eQ/Hs) + tan{eQ fBs')) = f dsL 

J —OO 


d is(p^—m^ 


-oo [dis 

. tan(eQ j: Bs) 2 

xe * {taTi{eQfBs) + tan{eQfBs')) (H.llO) 


Upon partial integration of Eqn. [H.110],i.e., 




{tan{eQ fBs) + tan{eQfBs'j) 


= — dse 

J —00 


dlS^ 


{tan{eQ fBs) + tan{eQ fBs')) 

2 -. tan(eq^Ba) 2 , , J 

— / dse ^11 “4/® ^ (tan(e(5/Hs) + tan(eQ/Hs')) (H.lll) 

J — 00 

It should be noted that, while deriving eqn. [H.lll] a total derivative term in has been thrown away, 
since the integrand vanishes at the boundary (virtue of the e prescription). Of the terms on the right 
hand side of Eqn. [H.lll], 

2 2 , ^ f .tan(eQj:Bs) 3 , , 'I 

dse'" PII * I (tan(eQ/Hs) + tan(e(5/Hs')) 

pOC) 

= (—1) / dse'^''''’"^p^sec^{eB) {tan{eQfBs) + taii{eQfBs')) 

J —00 


(H.112) 


and. 


^00 2 2 , .tan(e< 3 y: 8 a) 3 1 1 P 

dse"''ll {tan{eQfBs) + tan{eQfBs')) 

-nr, dlS 


/ OO 

e‘^^P'"hec^{eQfBs)ds. (H.113) 

-OO 


The expression for <l>(p, s) is the same as provided earlier. Therefore using Eqns.[ H.112] and [H.113 

/ OO poo 

ds — m^)(tan(eQ/Hs) + tan(e(5/Hs')) = / [p5_ sec^(eQ/H) (tan(e(5/Hs) 

-OO J —OO 

tan^(eQ/Hs')sec^(eQ/Hs) 2 


+tan(e(5/Hs'))+ieQ/Hsec^(e(5/H)] =/ 


(tan(eQ/Hs) + taii{eQ fBs')) 


(H.114) 
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Where in the last step of Eqn. [H.114], one has used the identity given by Eqn. [V.41]. It can now be 
easily verified, that on substituting Eqn. [H.114] in Eqn. [H. 109] the pieces proportional to k'j_ compensate 
each other. So the proof that, = 0 is complete. 
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